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PACS 42 . 50 . Pq - Cavity quantum electrodynamics; micromasers 

PACS 11.30.Er - Charge conjugation, parity, time reversal, and other discrete symmetries 
PACS 81.05.Xj - Metamaterials for chiral, bianisotropic and other complex media 

Abstract -Inspired by a recently observed asymmetry in the transmission of circularly polarized 
light through a metamaterial, we present a non-hermitian CPT-symmetric quantum model to de- 
scribe the interaction of the light fields in two resonant cavities coupled via a 2D-chiral mirror. We 
compute the time evolution of the light fields in this model, find two sets of operators compatible 
with the hamiltonian in a delocalized representation, discover the energies of the system and show 
that the transmission probability predicted by the model is indeed asymmetric. 
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Introduction. — Recently, a remarkable 4dB asym- 
metry was observed in the transmission of circularly polar- 
ized electromagnetic waves through an array of asymmet- 
ric split rings [1] . This implies that the total transmission 
probability for circularly polarized waves incident on the 
front face of the array is very different from the total trans- 
mission probability for circularly polarized waves incident 
&n the back face of the array. What would happen if this 
metamaterial array acted as a chiral semi-transparent mir- 
ror coupling two resonant cavities? We would expect the 
time evolution of the field in each cavity to be different 
from what is observed in the case of two cavities coupled 
through a regular mirror, but how different? 

In this paper, we propose a non-hermitian model that 
describes the quantum behaviour of the light field present 
in two resonant cavities coupled by such a 2D-chiral mir- 
ror. Firstly, we will resume some essential facts about 
optical cavities coupled by a standard, reciprocal interac- 
tion. Then, we will present our model for a non-reciprocal 
interaction between coupled cavities and discuss some of 
its consequences. 

Coupled cavities with reciprocal interaction. 

The description of one particular mode of the electromag- 
netic field in a cavity is equivalent to the description of 
a quantum oscillator. In this paper, we consider a very 



simple system consisting of two identical coupled quan- 
tum oscillators (A and B) with free ground-state energy 
Tiujq/2, described in the rotating- wave approximation by 
the hamiltonian 

H = H Q + H lth (1) 
with a free hamiltonian 



H = 1luj Q (a f a + 6 f 6 + l) 
and an interaction hamiltonian 



(2) 



(3) 



Each quantum oscillator represents the light field in one 
of the cavities. The coupling between the cavities is rep- 
resented by the interaction hamiltonian that describes the 
reciprocal interaction between two ordinarily coupled cav- 
ities. For definiteness, we choose g > 0. 

In terms of the quadrature operators (qa,Pa) for cavity 
A, and (qb,Pb) for cavity B, the Dirac operators are given 
by 

a = a W = ^T^A + ^a) b = b(t) = — ^=(?b + ipB) 
x> 2n x> 2n 



( a )E-mail: rsantosOfei.edu.br 



(4) 

a Ht) = —^(lA - ipA) b f = b\t) = -7=(qB - ipEs)- 
V2n v 2a 

(5) | 



p-1 



R. B. B. Santos 



We assume that the lowering operators a and b, and the 
raising operators and W all obey the usual equal-time 
canonical commutation relations 



[a, aJ] = 

We also assume that, at t 



1 = M f ] 
= o. 



a \n A , tib) = \/nA \ua - l,n B ) 
b \n A , n B ) = y/riB \n A , n B - 1) 
% \n A , n B ) = y/n A + 1 \n A + 1, n B ) 
% \n A , n B ) = \Jn B + 1 |ru, n B + 1) 



(0) 



(7) 
(8) 
(9) 
(10) 



where the ket | n A ,n B ) describes the state of the system 
in terms of the number of photons present in each cavity. 

The evolution of the light field is ruled by Heisenberg 
equations for the Dirac operators: 



[a,H] = -i(u a-gb) (11) 

[b,H] = -i(-ga + cj b) (12) 

[at,fl]=i(cjoo+- 5 6t) (13) 

[b\H] =i(- g a t +wo6 t ). (14) 



which may be solved by introducing lowering (pT and j3~) 
and raising (a + and /3 + ) operators for the two intercavity 
field modes by means of the canonical transformation 
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dt 
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~d7 





V2 



V2 



i-b 

7T 

xt - fct 



(15) 
(16) 



It should be remarked that the intercavity field modes 
are delocalized, consisting in a superposition of localized 
cavity modes from both cavities. It is a simple task to 
show that 



a 
a 1 



an e 



+ Acu a t 



/Jp-e-^ 4 
/3+e^*. 



where 



j a =uj - g top = Lu a + g 



(17) 
(18) 

(19) 



are the eigcnfrcquencies of the intercavity field modes a 
and P, respectively. Hence, the Dirac operators are given 
by 



a 



(cos(gt)a + i sin(gt)& ) e" iw °* (20) 

(ism(gt)a + cos{gt)b ) e" iw °* (21) 

(cos(gi)al - i sin(gt)&J)e ia;ot (22) 

(-isin(£/i)aj + cos(gt)&J)e iwot , (23) 



while the photon number operator for cavities A and £> 
are 

iV A = cos 2 (gt)a\a + sin 2 (gt)b\b 

+ icos(gt) s'm(gt)(alb — ao&J), (24) 

and 

7V B = sin 2 {gt)a\a Q + cos 2 (gt)blb 

+ icos{gt) sm(gt)(-alb + a bl), (25) 

respectively. 

Time evolution of the system is controlled by the unitary 
operator U = exp(— ii?t/7i). In the delocalized basis, 



U = exp(-iS„ ai „„t) \n a ,np) {n ai np\ 



(26) 



where E na>n „ = h(uj a n a -\-ujprif) In the localized 

basis, however, the time evolution operator is not diago- 
nal, leading to transitions between localized states. Since 
[Hq, i?i,h] = 0, the time evolution operator may be written 
as an operator product, 

IX = cxp(— iHot/h) cxp(— \Hifrt/h) 

« eM-iH t/h) (1 - iH^t/h) , (27) 

where the first order approximation suffices for the one- 
photon exchange between cavities. Apart from a phase 
factor contributed by exp(— iH^t/K), 

(n' a ,n' b \U\n a ,n b ) w 5™?'"? 

+ itg(yn a (n b + l)8<X-i + VK + iXiS). 

(28) 

characterizing a reciprocal interaction for the one-photon 
exchange probability between cavities A and £?. In the 
following sections, we present a model in which the two 
cavities are coupled by a non-reciprocal mirror and high- 
light some of its consequences. 

Description of the model. — The model we are 
proposing to describe the interaction of two cavities cou- 
pled by a non-reciprocal interaction is a variant of the 
coupled cavities model presented in section „ in which we 
change the interaction hamiltonian to 



I-h 



V.9AB 



SBACl 



(29) 



In this model, each quantum oscillator represents the light 
field in one of the cavities. The coupling between the cav- 
ities is represented by the interaction hamiltonian, which 
seeks to describe the transmission asymmetry aspects of 
the action of a 2Z?-chiral mirror. This interaction hamil- 
tonian is the simplest extension of the hermitian operator 
iJi.h that describes the reciprocal interaction between two 
ordinarily coupled cavities. We will restrict the present 
treatment to a closed system approach, disregarding any 
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interaction of the cavities with a reservoir. The interac- 
tion of the cavities with matter, and with a reservoir will 
be dealt with elsewhere. 

The hamiltonian H = Hq + Hi is a hermitian operator 
only if the condition <;ab = 3ba * s satisfied. In any other 
case, this hamiltonian will be non-hermitian. Particularly 
interesting is the case in which <7ab, <7ba G B,. In this case, 
even if <?ab 7^ <7ba, the hamiltonian H is a CRT-symmetric 
operator. In this non-reciprocal case, transmission proba- 
bilities are expected to be asymmetrical since <?ab 7^ 3ba- 

Parity (space inversion) CP and time reversal T are the 
basic operations of CPT symmetry. Under a parity trans- 
formation, q — > —q, and p — > — p while q — > q, p — > — p, 
and i — > — i under the anti- linear time reversal transforma- 
tion. Besides commuting with each other, the operators 
CP and T that realize parity and time reversal transforma- 
tions also satisfies CP 2 = 1 = T 2 . Since CP is unitary and 
T is anti-linear, the combined CPT operator is an antiuni- 
tary operator [2J. Under a CPT transformation, q — > —q, 
p — > p, and i — > — i. Hence, under a spacetime inversion 
operation, Dirac operators just reverse their signs, that is 
a, — > —a, — at. 

CPT-symmetric models have been receiving a great deal 
of attention since Bender and Boettcher [3] discussed a 
class of non-hermitian CPT-symmetric hamiltonian with 
real and positive eigenvalues. However, a CPT-symmetric 
hamiltonian alone is not enough to guarantee a real spec- 
trum. In fact, if the eigenvectors of the hamiltonian 
are not simultaneously eigenvectors of the CPT operator, 
then the lowest pair of eigenvalues dissolves into com- 
plex conjugates, the spectrum ceases to be real and the 
CPT symmetry of the hamiltonian is said to be broken [3] . 
This spontaneous CPT symmetry breaking, and the excep- 
tional points [SJH] or the non-hermitian degeneracies [7] 
at which the symmetry breaking occurs have been the 
subject of intense research in recent years. Even optical 
valves and cloaking devices based on CPT symmetry break- 
ing were recently proposed [H[5] and analysed [TU] . A few 
years after the pioneer work by Bender and Boettcher, 
Mostafazadch [TT] showed that exact CPT-symmetry is 
equivalent to hermiticity with respect to a suitably de- 
fined inner product. Recent reviews [12l[T3] provide a de- 
tailed exposition of the status of the theory. Experimental 
realizations of systems with CPT-symmetry appeared much 
more recently, initially in the optics domain pBHTS] follow- 
ing a theoretical analysis of coupled optical CPT-symmetric 
structures involving a balanced gain- loss profile [17] . More 
recently, an active RLC circuit realization was demon- 
strated [IS]- Despite the impressive results achieved in 
these experiments, they all share the common trait of 
appearing to be strictly classical systems. However, a 
theoretical analysis [19ll20] based on the scattering for- 
mulation [2TH23] of the input-output formalism [24, 25 
in quantum optics showed that quantum noise, losses 
through leakage, and instability are critical to compre- 
hend the observed behaviour of these systems in which 
non-hermiticity arises from coupling lossy and active el- 



ements, that is, absorbing and amplifying regions, via a 
reciprocal coupling. In our model, however, the source 
of non-hermiticity is the presence of a non-reciprocal cou- 
pling between two otherwise ordinary resonant cavities. 
Non-reciprocal Bragg scattering from a photonic lattice 
with CPT symmetry was predicted |26j . owing to a break- 
down of invariance of the diffracted beam intensities under 
crystal inversion [2"Tll2"8] . Still relating to CPT symmetry in 
optical cavities, Longhi [29] built on the ideia of a coher- 
ent perfect absorber [30] and used a classical scattering 
formalism to show that a CPT-symmetric optical medium 
in a resonant cavity may behave both as a laser oscillator, 
and as a coherent perfect absorber. Another interesting 
development is the search for optical realisations of rela- 
tivistic non-hermitian systems [31] , We believe that the 
system analysed in this paper can be realised combining 
standard techniques from cavity quantum electrodynamics 
and metamaterials research. Therefore, it could be one of 
the simplest experimental realizations of a CPT-symmetric 
quantum system. 

In the following sections, we will explore the conse- 
quences of this non-hermitian hamiltonian, owing to the 
non-reciprocal nature of the interaction in the case |<?ab| 7^ 
I <7ba I ■ This case is clearly impossible in a hermitian model. 

Equations of motion. — The dynamics of the system 
is ruled by Heiscnbcrg equations for the evolution of the 
Dirac operators: 



da 
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dt 


~ m 


db 
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at 


~ m 


dot 
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dt 


~ m 


dtf 


1 


~dt 





[a, H] = -i {u> Q a - g B Ab) 
[b, H] = -i (-ffABa + uj b) 



(30) 
(31) 
(32) 



It should be noticed that 



dot /da x 1 

~dT ^ \ dt 

d&t /db\ 1 

~dT^{di 



(34) 
(35) 



that is, lowering and raising operators behave as truly in- 
dependent dynamical variables in this model. This con- 
trasts with the hermitian formulation, in which the two 
Dirac operators associated with each quantum oscillator 
are just hermitian conjugates of each other. This result 
is not really surprising since CPT symmetry is weaker than 
hermiticity. 

The two sets of coupled operator differential equa- 
tions ([3"D"]) and (|3"3"]) may be solved by the canonical trans- 
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formation 



o ' 



V 2 5ba 




AB 



AB 



(37) 



These new operators are the lowering (a~ and (3~) and 
the raising (a + and /3 + ) operators for the two intercav- 
ity field modes. It should be remarked, however, that the 
intcrcavity field modes are delocalized, consisting of a su- 
perposition of localized cavity modes from both cavities. 
Nonlocal hermitian interactions may, in some cases, be de- 
scribed more simply by a non- hermitian hamiltonian . 
From now on, we will refer to the description of the system 
given by the cavity field mode operators a, a* , b, and W as 
the localized representation while the expression delocal- 
ized representation will be associated to the description of 
the system given by the intercavity field mode operators 
a~ , a + , (3~ and (5 + . 

Time evolution in the delocalized representation. 

— The intercavity field mode operators satisfy the ex- 
pected equal-time canonical commutation relations and 
their evolution is ruled by the uncoupled differential equa- 
tions 



dt 

d + • -i 
dt 



iw (3 /J + , 



(38) 
(39) 



where 



uj a = u - \/.9ab.9ba up = ujq + VSABffBA (40) 

are the eigenfrequencies of the intercavity field modes 

a and /?, respectively. Equations (|3"5j) - (|3"9"]) are solved 
straightforwardly: 



n 



= a e 



+ ioj a t 



a 6 e 



r 



/?o + e 



(41) 
(42) 



In the delocalized representation the hamiltonian is 
quite simple: 



H = H (cv a N a + upN? + w ) 



(43) 



where the excitation number operators for the intercavity 
field modes are given by 



a o a o 



N a = a+a 



(44) 
(45) 



As expected, the number of excitations in each intercavity 
field mode docs not depend on time. 

Since the excitation number operators for the intercav- 
ity field modes are compatible observables, \n a ,np) is an 
eigenstate of the hamiltonian with energy 



h (uj a n a + LUprip + uj ) 



(46) 



In the regime of weak to moderate intercavity coupling 
(\/.9ab.9ba < ^o), energy is always positive. However, in 
the strong intercavity coupling regime (-\/<?ab9ba > wo), 
En ai np m ay be negative, signalling a breakdown of the 
rotating-wave approximation. Beside the set {N a , N"}, 
another complete set of operators for the system is formed 
by the total excitation number operator N = N a +N@ and 
by the excitation imbalance operator A Q( 3 = N a — N 13 . In 
terms of these operators, the hamiltonian may be cast in 
the form 



H = hajQ (N + 1) - fiV5AB5BAA Q/3 . 



(47) 



We could not find any evidence of a broken yT-symmetry 
phase in the spectrum of the eigenvalues of H . In fact, 
it is not clear whether or not this system is required to 
present a broken [PT-symmetry phase as it does not seem 
to fit in any of the known universality classes discussed by 
Schomerus [52"] . 

Although the delocalized representation based on intcr- 
cavity field mode operators is very useful in the theoretical 
description of the coupled system, the delocalized observ- 
ables required may not be very easy to measure compared 
to their localized counterparts iV A = a^a and N B = b^b 
which represent the number of photons in each cavity. 

In terms of the localized cavity operators, the total exci- 
tation number and the excitation imbalance operators are 
simply 



N = aja + b\b a 

A a/ 3 = 



3BA ti, 

%bo 

9AB 



SAB ,t 

a %. 

9BA 



(48) 
(49) 



In the next section, we will investigate the time evolution 
of the operators of the localized representation. 

Time evolution in the localized representation. 

The cavities lowering and raising operators are determined 
from equations (|4"T|) ~ P2]) 



cos(^3AB3BAi)ao + l \J^^ sm (V.9ABffBAi) fo o j e 

(50) 



( ~~ sin(^/3ABffBA*)ao + cos(^.g A B.gBA£)fro ) e " 



,t _ 



cos(^3AB5BAt)aJ - L / sin(^/g A Bfl , BA*)&o ) e 



(51) 



t | Ju t 



9BA 



(52) 



= (^-y— sm(y/ g AB g B At)al + cos^abSba^oJ e 1 " ' 

(53) 

As a consistency check, we should point out that the equal- 
time canonical commutation relation © still holds at any 
time in this non-hermitian model with non-reciprocal in- 
teraction. 
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Within the localized representation, the photon number 
operators iV A = a^a and iV B = rtb for cavities A and B, 
respectively, are incompatible with the hamiltonian, and, 
therefore, varies in time, as expected. It is simple to show 
that 

N A = cos 2 (y.9AB5BAi)aS a o + sm 2 (^/g AB g BA t)blb 
+ . cos( % /5AB.9BAi) sin(^/g AB g BA t)x 

V5ABffBA 

(gBAalbo - g A Baobl), (54) 

and that 

N B = sin 2 (^9Aa3BAt)aJa + cos 2 (^5 A B5BAi)fy^o 
+ . cos(^.g A B5BAi) sin( v /pABfl'BA*)x 

V5AB.9BA 

(-.9BAa^o + gABaobl). (55) 

Besides, the average photon number in the cavities be- 
haves as expected. For cavity A, we have 

(n A ,n B \ N A \n A ,n B ) = n A cos 2 (^/g AB g BA t) 

+ n B sin 2 (^/g A B.9BAi) (56) 

while, for cavity B, 

(n A ,n B \ N B \n A ,n B ) = n A sin 2 (y/g AB g BA t) 

+ n B cos 2 {^/g AB g BA t). (57) 

In fact, since the system is not coupled to a reservoir, 

N = N A + N B =N A + N B = ojoo + blb , (58) 

that is, at any time, the number of photons in the coupled 
cavities is always equal to the initial number of photons in 
the system, as expected. The same conclusion is achieved 
from an analysis of equation (|4"5)) for the total excitation 
number operator. 

Computation of probability amplitudes. In the 

localized representation, the action of the lowering opera- 
tor a for cavity A on a localized number state Itiaj^b) is 
given by 

a|n A ,n B )= I \/n A cos(^/g AB g BA t) |n A - l,n B ) 



- H / — y/n£ sin( V<?ab <?ba*) |n A ,n B - 1) c^ ' (59) 

V 9AB J 

whereas the action of the raising operator at for cavity A 
is 



a) \n A ,n B ) = \^/n A + 1 cos(^.g A B.9BAi) \n A + l,n B ) 
+ uf^Vm + lsm(V9AB9BAt) |n A ,n B + 1) )c iu "' t 

V .9BA / 

(60) 



which are both quantum superpositions in the localized 
representation. Similar equations hold for b, and for rt. 
In fact, it is worth pointing that 

Hi |nA,n B ) = -%ab \/n A {n B + 1) \n A - l,n B + 1) 
- hg BA \J (n A + l)n B \n A + 1, ?i B - 1) . 

(61) 

As a consistency check, we note that the vacuum state 
| vac) = 1 0,0) is preserved by the interaction hamiltonian. 
Another interesting situation arise when the initial state 
is a pure cavity number state characterized by the same 
number of photons in both cavity. In this case, the inter- 
action hamiltonian maps the state \n, n) into an entangled 
state of the |n + 1, n — 1) and \n — 1, n + 1) states. 

Since the free hamiltonian Hq commutes with the inter- 
action hamiltonian Hi, we may use the Baker-Campbell- 
Hausdorff formula to obtain 

It = exp(-iHt/h) = exp(-iH t/h) exp(-iHxt/h). (62) 

Unfortunately, 

[ab\a%} = N B - N A (63) 

is not a scalar. Since 

[art, [ab\ at b]] = -2art (64) 
[a%, [ab\a^b]] = 2atb, (65) 

we would have to consider higher order terms in the Baker- 
Campbcll-Hausdorff formula, which is cumbersome. We 
turn for the small-time expansion of cxp(—iHit/h) « 
1 — iHit/h around the initial time t = 0, instead. This 
approximation suffices in order to draw a physically mean- 
ingful picture of the interaction process. 

Apart from a phase, the matrix elements of the evolu- 
tion operator U are 

(n' A ,n' B \U\n Al n B ) « 5„ A ,n A £n B ,n B 

- it(g AB ^n A (n B + 1)<5„ a ,„ a -iC>„ b ,„ b + 



+ .obaV ("a + l)n B (5„ A ,„ A +i<5„^,„ B _i 

(66) 

Close inspection of Eq. (|6"6"|) reveals that the probabil- 
ity amplitude for a photon to be transmitted from cavity 
A to cavity B is proportional to g AB \/n A (n B + 1) while 
the probability amplitude for a photon to be transmit- 
ted back from cavity B to cavity A is proportional to 
ffBA\/ (i^a + 1)?t-b- Since these amplitudes are, in general, 
not the same, we may say that the quantum system de- 
scribed in this paper is asymmetrical regarding the ex- 
change of photons. This asymmetry is a consequence of 
CRT-symmetry in a quantum system. It must be remarked 
that this kind of asymmetry is impossible in a closed sys- 
tem described by a hermitian hamiltonian. 
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Discussion. — We proposed a non-hermitian 77- 
symmetric model for describing the light field of two cavi- 
ties interacting through a 2-D-chiral mirror. Analysing the 
time evolution of the system, we showed that the probabil- 
ity amplitude for the exchange of a photon from one cavity 
to another through the 2_D-chiral mirror differs from the 
probability amplitude for the reverse process. This re- 
sult suggests a new possible application of non-hermitian 
quantum theories, specially CPT-symmetric quantum me- 
chanics, as an effective theory dealing with some quantum 
aspects of light propagation through optical rn.etam.ate- 
rials. Besides, we must point out the fact that, in this 
model, the Dirac operators behave as truly independent 
degrees of freedom, and not just the hermitian conjugates 
of each other. 

* * * 

We acknowledge an anonymous referee for comments 
that improved the paper and Fr. Saboia de Medeiros Ig- 
natian Educational Foundation for continuous support. 
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